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Updated Estimate of the Muon Magnetic Moment 
Using Revised Results from e + e _ Annihilation 
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A new evaluation of the hadronic vacuum polarization contribution to the muon magnetic moment is presented. 
We take into account the reanalysis of the low-energy e + e~ annihilation cross section into hadrons by the CMD-2 
Collaboration. The agreement between e + e~ and r spectral functions in the tt-k channel is found to be much 
improved. Nevertheless, significant discrepancies remain in the center-of-mass energy range between 0.85 and 
1.0 GeV/c 2 .The deviations from the measurement at BNL are found to be (22.1 ± 7.2 ± 3.5 ± 8.0) 10~ 10 (1.9 cr) 
and (7.4 ± 5.8 ± 3.5 ± 8.0) 10 -10 (0.7 a) for the e + e~- and r-based estimates, respectively, where the second error 
is from the nonhadronic contributions and the third one from the BNL measurement. Taking into account the 
p~ — p° mass splitting determined from the measured spectral functions increases the r-based estimate and leads 
to a worse discrepancy between the two estimates. 



1. Introduction 

Hadronic vacuum polarization in the photon 
propagator plays an important role in the pre- 
cision tests of the Standard Model. This is the 
case for the muon anomalous magnetic moment 
a /j = — 2)/2 where the hadronic vacuum po- 
larization component, computed from experimen- 
tally determined spectral functions, is the leading 
contributor to the uncertainty of the theoretical 
prediction. 

Spectral functions are directly obtained from 
the cross sections for e + e~ annihilation into 
hadrons. The accuracy of the calculations has 
therefore followed the progress in the quality of 
the corresponding data pQ. Because the latter 
was not always suitable, it was deemed neces- 
sary to resort to other sources of information. 
One such possibility was the use of the vector 
spectral functions [2] derived from the study of 
hadronic r decays J5J for the energy range less 
than m T ~ 1.8 GeV. Also, it was demonstrated 
that perturbative QCD could be applied to en- 
ergy scales as low as 1-2 GeV @], thus offering 
a way to replace poor e + e~ data in some energy 



regions by a reliable and precise theoretical pre- 
scription |5I6I7I8I9] . 

A complete analysis including all available ex- 
perimental data was presented in Ref. jlU) . tak- 
ing advantage of the new precise results in the 
7T7T channel from the CMD-2 experiment JT] 
and from the ALEPH analysis of r decays |12) . 
and benefiting from a more complete treatment 
of isospin- breaking corrections |13I14) . In addi- 
tion to these major updates, the contributions 
of the many exclusive channels up to 2 GeV 
center-of-mass energy were completely revisited. 
It was found that the e + e~ and the isospin- 
breaking corrected r spectral functions were not 
consistent within their respective uncertainties, 
thus leading to inconsistent predictions for the 
lowest-order hadronic contribution to the muon 
magnetic anomaly. The leading contribution to 
the discrepancy originated in the 7T7T channel 
with a difference of (— 21.2 ± 6.4 oxp ± 2.4 rac i ± 
2-6su(2) (±7.3 to tai)) 10~ 10 . The estimate based 
on e + e~ data agreed with another analysis us- 
ing the same input data |15| . When compared to 
the world average of the muon magnetic anomaly 



2 



measurements, 

a M = (11 659 203 ±8) 1(T 10 , (1) 

which is dominated by the 2002 BNL result using 
positive muons {H^, the respective e + e~-based 
and r-based predictions disagreed at the 3.0 and 
0.9 a level, respectively, when adding experimen- 
tal and theoretical errors in quadrature. 

Our analysis had to be updated ^7] since the 
CMD-2 Collaboration at Novosibirsk discovered 
that part of the radiative treatment was incor- 
rectly applied to the data and produced a com- 
plete reanalysis As the CMD-2 data domi- 
nate the e + e _ -based prediction, the changes pro- 
duce a significant effect in the final result. No 
significant change occurred for the r-based pre- 
diction. The only relevant fact is a new result [17?] 
for the branching ratio of the r~ — ► v T h~ir° mode 
(h~~ stands for a charged pion or kaon). 

2. Muon Magnetic Anomaly 

It is convenient to separate the Standard Model 
(SM) prediction for the anomalous magnetic mo- 
ment of the muon into its different contributions, 

*T = a^+a^ + a^\ (2) 

with 

had had.LO , had, HO , had.LBL /o\ 

a p = % ' + a p + a p , (3) 

and where a^ ED = (11 658 470.6 ± 0.3) 10~ 10 is 
the pure electromagnetic contribution (see |20I21) 
and references therein : ), a h ad > LO is the lowest- 
order contribution from hadronic vacuum polar- 
ization, a h ad . H o _ (_io.O±0.6) 10" 10 is the cor- 
responding higher-order part |24I2| . and a™ cak = 
(15.4 ± 0.1 ± 0.2) 10" 10 , where the first error 
is the hadronic uncertainty and the second is 
due to the Higgs mass range, accounts for cor- 
rections due to exchange of the weakly interact- 
ing bosons up to two loops [35] ■ For the light- 
by-light (LBL) scattering part we add the val- 
ues for the pion-pole contribution 26 27 28 and 

1 Some adjustment was recently made concerning the 
fourth-order contribution from the leptonic light-by-light 
scattering, mostly affecting the QED prediction for a e and 
through it the value of a 1221231 . The resulting change in 
a^ ED is within the quoted uncertainty of 0.3 10~ 10 and 
has not been included in the present analysis. 



the other terms |27I28| to obtain a Jf d ' LBL = 
(8.6 ±3.5) 10- 10 . 

Owing to the analyticity of the vacuum po- 
larization correlator, the contribution of the 
hadronic vacuum polarization to can be cal- 
culated via the dispersion integral |29"] 

oo 

«r LO = ^ Jd S ^R( S ), ( 4) 

where K(s) is a well-known QED kernel. In 
Eq. (gj, R(s) = i? (0) (s) denotes the ratio of the 
'bare' cross section for e + e~ annihilation into 
hadrons to the pointlike muon-pair cross section. 
The 'bare' cross section is defined as the mea- 
sured cross section, corrected for initial-state ra- 
diation, electron-vertex loop contributions and 
vacuum polarization effects in the photon prop- 
agator (note that photon radiation in the final 
state (FSR) is included in the 'bare' cross sec- 
tion). The reason for using the 'bare' (i.e. low- 
est order) cross section is that a full treatment 
of higher orders is anyhow needed at the level 
of a M , so that the use of 'dressed' cross sections 
would entail the risk of double-counting some of 
the higher-order contributions. 

The function K(s) decreases monotonically 
with increasing s. It gives a strong weight to the 
low energy part of the integral QJ. About 91% of 
the total contribution to a^ad.LO j g accum ulated 
at center-of-mass energies y/s below 1.8 GeV and 
73% of a^ ad ' LO is covered by the two-pion final 
state which is dominated by the /?(770) resonance. 

3. Changes to the Input Data 

The CMD-2 data, published in 2002 for the 
irir channel |11| . have been completely reana- 
lyzed PB] following the discovery of an incor- 
rect implementation of radiative corrections in 
the analysis program. Overall, the pion-pair cross 
section increased by 2.1% to 3.8% in the mea- 
sured energy range (cf. Fig. QJ, well above the 
previously quoted total systematic uncertainty of 
0.6%. Specifically, the leptonic vacuum polariza- 
tion contribution in the t-channel had been inad- 
vertently left out in the calculation of the Bhabha 
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Figure 1. Relative change in the e + e~ — > tt + h~ 
cross section of the revised CMD-2 analysis JTffl 
with respect to the one previously published lllf . 



cross section. This effect produced a bias in the 
luminosity determination, varying from 2.2% to 
2.7% in the 0.60-0.95 GeV energy range. The 
problem consequently affected the measured cross 
sections for all hadronic channels. Another prob- 
lem was found in the radiative corrections for the 
muon-pair process, ranging from 1.2% to 1.4% 
in the same region. A more refined treatment 
of hadronic vacuum polarization was performed, 
with changes not exceeding 0.2% for most data 
points. The effects in the Bhabha- and muon- 
pair channels also affected the event separation 
and the measured ratio of pion pairs to electron 
and muon pairs changed by typically 0.7%. 

The correction of the bias in the luminosity de- 
termination increases all hadronic cross sections 
published by CMD-2. The changes are 2.4% and 
2.7% on the u> and resonance cross sections, re- 
spectively. 

New published data by SND on the u> reso- 
nance 1201 and the 2-k + 2tt~ as well as tt + tt~2tt° 
modes jHJ (unchanged cross sections for the lat- 
ter two, but reduced systematics with respect to 
previous publications) have been included in this 
update. 

A detailed discussion of radiative corrections, 
in particular the effect of final-state radiation 
by the charged hadrons was given in Ref. |10) . 



Also given therein is a compilation of all input 
data (with references) used to calculate the inte- 
gral (H}. 

4. Comparison of e + e and r Spectral 
Functions 

The new e + e~ and the isospin-breaking cor- 
rected t spectral functions can be directly com- 
pared for the 7T7T final state. The r spectral 
function is obtained by averaging ALEPH [3], 
CLEO [12 and OPAL 33; results PHI . The e+e" 
data are plotted as a point-by-point ratio to the r 
spectral function in Fig.0 and enlarged in Fig. El 
to better emphasize the region of the p peak. The 
central bands in Figs. |21 and give the quadratic 
sum of the statistical and systematic errors of the 
r spectral function obtained by combining all r 
data. The e + e~ data have moved closer to the r 
results: they are now consistent below and around 
the peak, while, albeit reduced, the discrepancy 
persists for energies larger than 0.85 GeV. 

A convenient way to assess the compatibility 
between e + e~ and r spectral functions proceeds 
with the evaluation of r decay fractions using the 
relevant e + e~ spectral functions as input. All the 
isospin-breaking corrections detailed in Ref. [TH] 
are included. This procedure provides a quan- 
titative comparison using a single number. The 
weighting of the spectral function is however dif- 
ferent from the vacuum polarization kernels. Us- 
ing the branching fraction B(r~ — ► v T e~ u e ) = 
(17.810 ± 0.039)%, obtained assuming leptonic 
universality in the charged weak current |12j , the 
result for the nn channel is 

B^c = (24.52±0.26 cxp ±0.11 rad ±0.12 su(2) )% , (5) 

where the errors quoted are split into uncertain- 
ties from the experimental input (the e + e~ an- 
nihilation cross sections) and the numerical inte- 
gration procedure, the missing radiative correc- 
tions applied to the relevant e + e~ data, and the 
isospin-breaking corrections when relating r and 
e + e~ spectral functions. Even though the cor- 
rections to the CMD-2 results have reduced the 
discrepancy between JSJ and the world average of 
the direct B(t~ — > zy T 7r~7r°) measurements from 
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24.52 + 0.32 



CLEO 

25.42 + 0.12 + 0.42 



OPAL 

25.44 + 0.17 + 0.29 



L3 
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ALEPH preliminary 
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t Average 
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Figure 2. Relative comparison of the 7r + 7r~ spec- 
tral functions from e + e~ and isospin-breaking 
corrected t data, expressed as a ratio to the r 
spectral function. The band shows the uncertainty 
on the latter. The e + e~ data are from CMD- 
2 \WI, CMD, OLYA and DM1 (references quoted 
in \T7j ). 



Figure 4. The measured branching ratios for 
t~ — > v T ir~ir° compared to the prediction from 
the e + e~ — > tt + it~ spectral function applying the 
isospin-breaking correction factors discussed in 
Ref. The measured branching ratios are from 

ALEPH CLEO '3^ and OPAL ]EB- The 
L3 and OPAL results are obtained from their hir° 
branching ratio, reduced by the small Kir contri- 
bution measured by ALEPH JSfijl and CLEO WH - 



• CMD-2 
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o CMD 
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■ DM1 
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4.6 to 2.9 standard deviations (adding all errors in 
quadrature), the remaining difference of (— 0.94± 
0.10 T ±0.26 O e±0.11 rad ±0.12 su(2) (±0.32 total ))%is 
still problematic. Since the disagreement between 
e + e~ and r spectral functions is more pronounced 
at energies above 850 MeV, we expect a smaller 
discrepancy in the calculation of a^ ad,LO because 
of the steeply falling function K(s). More infor- 
mation on the comparison is displayed in Fig. 0] 
where it is clear that ALEPH, CLEO, L3 and 
OPAL all separately, but with different signifi- 
cance, disagree with the e + e~-based CVC result. 



Figure 3. Relative comparison in the p region 
of the 7r + 7r _ spectral functions from e + e~ and 
isospin-breaking corrected r data, expressed as a 
ratio to the r spectral function. The band shows 
the uncertainty on the latter. The references for 
e + e~ data are given in Fig. 



5. Results 

The integration procedure and the specific con- 
tributions - near tttt threshold, the tu and <fi reso- 
nances, the narrow quarkonia and the high energy 
QCD prediction - are treated as in our previous 
analysis JUj. The contributions from the differ- 
ent processes in their indicated energy ranges are 
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Not including w and c/i resonances (see Ref. [To|1. 
2 Using isospin relations (see Ref. HUl V 
3 -.e+e — data are used above 1.6 GeV (see Ref. IIUI I. 



Summary of the a^ ad ' LO contributions from e + e~ annihilation and r decays. The uncertainties on the 
vacuum polarization and FSR corrections are given as second errors in the individual e + e~ contribu- 
tions, while those from isospin breaking are similarly given for the r contributions. These 'theoretical' 
uncertainties are correlated among all channels, except in the case of isospin breaking which shows little 
correlation between the 2n and 4V channels. The errors given for the sums in the last line are from the 
experiment, the missing radiative corrections in e + e~ and, in addition for t , SU(2) breaking. 
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listed in Table ^ Wherever relevant, the two 
e + e~- and r-based evaluations are given. The 
discrepancies discussed above are now expressed 
directly in terms of a^ ad ' LO , giving smaller esti- 
mates for the e + e~-based data set by (—11.9 ± 
6.4 cxp ± 2.4 rad ± 2.6 SU(2) (±7.3 t otai)) KT 10 for 
the 7T7T channel and (—2.8 ± 2.6 cxp ± 0.3 ra d ± 
1.0 SU(2 ) (±2.9 to tai)) 10~ 10 for the sum of the 4?r 
channels. The total discrepancy (— 14.7±6.9 exp ± 
2.7 rad ±2.8 su(2) (±7.9 to tai)) 10~ 10 amounts to 1.9 
standard deviations. The difference could now be 
considered to be acceptable, however the system- 
atic difference between the e + e~ and r 7T7t spec- 
tral functions at high energies precludes one from 
performing a straightforward combination of the 
two evaluations. 

5.1. Results for a /j, 

The results for the lowest order hadronic con- 
tribution are (xlO~ 10 ) 

<? e d e L ° = 696.3 ± 6.2 cxp ± 3.6 rad (6) 
fl had,LO = rii.o ±5.0 exp ±0.8 rad ±2.8 su(2 ^7) 

Adding to these the QED, higher-order hadronic, 
light-by-light scattering and weak contributions 
as given in Section [3 we obtain the SM predic- 
tions (x 10~ 10 , with an additional common uncer- 
tainty of ±3.5 LB l ± 0.4 QED+EW ) 

a™ e = 11 659 180.9 ±7.2 had xo 

a™ = 11 659 195.6 ±5.8 h ad,LO 

These values can be compared to the present mea- 
surement Adding experimental and theoret- 
ical errors in quadrature, the differences between 
measured and computed values are found to be 
(x 10 -10 , with an additional common uncertainty 
of ±3.5 o ther±8.0 oxp from contributions other than 
hadronic vacuum polarization and the BNL g-2 
experimental error) 

a c * p ~al M ee = 22.1±7.2 had , LO (] 
< p ~a™ = 7.4 ± 5.8 had , LO [ ) 

where the error quoted is specific to each ap- 
proach. The differences (JBJ) correspond to 1.9 and 
0.7 standard deviations, respectively. A graphi- 
cal comparison of the results © with the exper- 
imental value is given in Fig. [S] Also shown are 
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Figure 5. Comparison of the results @) with 
the BNL measurement }lftf - Also given are our 
estimates f-?l£f obtained before the CMD-2 data 
were available. For completeness, we show as tri- 
angles with dotted error bars the e + e~ -based re- 
sults 11 (Ml 51/ derived with the previously published 
CMD-2 data HI}/ . 



our estimates |1I9| . obtained before the CMD-2 
and the new r data were available (see discus- 
sion below), and the e + e _ -based evaluations of 
Rcfs. [10115) . obtained with the previously pub- 
lished, uncorrected CMD-2 data [11| . 

6. Discussion 

Although the new corrected CMD-2 ir + ir~ re- 
sults are now consistent with r data for the energy 
region below 850 MeV, the remaining discrep- 
ancy for larger energies is unexplained at present. 
Hence, one could question the validity of either 
e + e _ data with their large radiative corrections, r 
data, or the isospin-breaking corrections applied 
to t data. We shall briefly discuss these points 
below. 

• The CMD-2 experiment is still the only one 
claiming systematic accuracies well below 
1%. It is thus difficult to confront their 
data with results from other experiments. 
Whereas the measurements from OLYA are 
systematically lower than the new CMD-2 
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results in the peak region, there is a trend of 
agreement above, as seen in Fig. 0] This be- 
haviour appears to be confirmed by prelim- 
inary data from the KLOE experiment at 
Frascati using the radiative return method 
from the <f> resonance |38| . 

• The most precise results on the r titt 
spectral function come from the ALEPH 
and CLEO experiments, operating in com- 
pletely different physical environments. On 
the one hand, the main uncertainty in 
CLEO originates from the knowledge of the 
relatively low selection efficiency, a conse- 
quence of the large non-r hadronic back- 
ground, while the mass spectrum is mea- 
sured with little distortion and good resolu- 
tion. On the other hand, ALEPH has both 
large efficiency and small background, the 
main uncertainty coming from the ir° recon- 
struction close to the charged pion, neces- 
sitating to unfold the measured spectrum 
from detector resolution and acceptance ef- 
fects. A comparison of the r spectral func- 
tions from ALEPH, CLEO and OPAL is 
given in Fig. EJ and in Fig. for the p 
peak region. Agreement is observed within 
quoted errors, in particular in the high mass 
region, although CLEO results are a bit 
closer to e + e~ data there. Overall, the r 
data appear to be consistent. 

• The last point concerns isospin corrections 
applied to the r spectral functions. The 
basic components entering SU(2) breaking 
are well identified. The long-distance radia- 
tive corrections and the quantitative effect 
of loops have been addressed by the anal- 
ysis of Ref. 23] showing that the effects 
are small. The overall effect of the isospin- 
breaking corrections (including FSR) ap- 
plied to the t 7T7T data, expressed in relative 
terms, is (—1.8 ± 0.5)%. Its largest con- 
tribution (—2.3%) stems from the uncon- 
troversial short-distance electroweak correc- 
tion PHj. One could question the validity 
of the chiral model used. The authors of 
Ref. 23] argue that the corrections are in- 
sensitive to the details of their model and 
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Figure 6. Relative comparison of the 7r + 7r~ spec- 
tral functions extracted from r data from different 
experiments, expressed as a ratio to the average t 
spectral function. 
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Figure 7. Relative comparison in the p region 
of the Tr + n~ spectral functions extracted from r 
data, expressed as a ratio to the average r spec- 
tral function. 
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essentially depend only on the shape of the 
pion form factor. As the latter is known 
from experiment to adequate accuracy, it 
seems difficult to find room for a ~ 10% 
effect as observed experimentally. Never- 
theless, considering the situation regarding 
the first two experimental points, it would 
seem worthwhile to invest more theoretical 
work into the problem of isospin breaking. 

The particular point of the p~ — p° mass 
splitting deserves some further discussion. 
This possibility was considered from the 
beginning when it was proposed to use r 
spectral functions to compute vacuum po- 
larization integrals [2]. However, at that 
time, the experimental investigation led to a 
splitting consistent with within 1.1 MeV, 
supported by theoretical investigations |3U] 
which indicated a value less than 0.7 MeV. 
Consequentlyin further analyses, we as- 
sumed m„o = m p - to hold within an un- 
certainty of 1 MeV. A re-analysis of this 
question with the current more precise data 
on e + er and r spectral functions leads 
to the conclusion that the mass splitting 
favoured by the data is m p - —m„o = (2.3 ± 

0. 8 MeV. A similar conclusion was more 
recently reached 021 • Unfortunately, it is 
not yet clear whether such a result must 
be taken as definitive evidence for a mass 
splitting and if the corresponding correction 
must be applied to the r spectral function. 
The significance of the result is still lim- 
ited, but the most worrisome problem is the 
fact that the e + e~ and the mass-corrected 
t spectral functions still disagree in magni- 
tude by 3.3% Therefore, while the con- 
sideration of the mass splitting improves the 
line shape comparison, it leaves us with a 
major normalization discrepancy. In other 
words, the p~ mass correction spreads the 
large local difference observed in the 0.85- 

1. GeV range almost uniformly across the 
mass spectrum. We thus disagree with the 
conclusions reached in Ref. 021 that taking 
into account the p mass difference solves 
the current discrepancy. If we take the 



point of view that the remaining difference 
is basically a normalization problem (either 
from the data or the isospin breaking cor- 
rections), then correcting for the apparent 
mass difference increases the r-based esti- 
mate of a had,LO ky 5 4 10~ 10 , bringing it 
even closer to the BNL result, but further 
away from the e + e~-based estimate for a 
total of 2.5 a. More experimental investi- 
gation with CMD-2, KLOE and BABAR is 
needed to consolidate the e + e~ picture. 



7. Conclusions 

An update of our analysis of the lowest-order 
hadronic vacuum polarization contribution to the 
muon anomalous magnetic moment has been per- 
formed following a reevaluation by the CMD-2 
Collaboration of their e + e~ annihilation cross 
sections. Part of the previous discrepancy be- 
tween the e + e~ and r irir spectral functions has 
now disappeared so that the corresponding eval- 
uations of the lowest-order hadronic polarization 
contribution to the muon magnetic anomaly are 
closer. However, incompatible measurements re- 
main between 0.8 and 1 GeV so that we do not 
proceed with an average of the two evaluations. 
The e + e _ - and r-based predictions are respec- 
tively 1.9 and 0.7 standard deviations below the 
direct measurement from the g-2 Collaboration at 
BNL. Considering the tttt discrepancy from the 
point of view of the p~ — p° mass splitting turns 
out to increase the difference between the two es- 
timates. The forthcoming results from radiative 
return with KLOE jHE] and BABAR gSj will be 
important to sort out possibly remaining prob- 
lems in the tttt and Air spectral functions. 
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